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Abstract 

As a first step towards a strong coupling expansion of 
Yang-Mills theory, the SU{2) Yang-Mills quantum me- 
chanics of spatially constant gauge fields is investigated 
in the symmetric gauge, with the six physical fields rep- 
resented in terms of a positive definite symmetric 3x3 
matrix S. Representing the eigenvalues of S in terms of 
elementary symmetric polynomials, the eigenstates of the 
corresponding harmonic oscillator problem can be calcu- 
lated analytically and used as orthonormal basis of trial 
states for a variational calculation of the Yang-Mills quan- 
tum mechanics. In this way high precision results are ob- 
tained in a very effective way for the lowest eigenstates 
in the spin-0 sector as well as for higher spin. Further- 
more I find, that practically all excitation energy of the 
eigenstates, independently of whether it is a vibrational 
or a rotational excitation, leads to an increase of the ex- 
pectation value of the largest eigenvalue (fo), whereas the 
expectation values of the other two eigenvalues, (</>i) and 
(</>2), and also the component (B3) — g{4>\4>2) of the mag- 
netic field, remain at their vacuum values. 

1 Introduction 

The so-called Yang-Mills mechanics originates from Yang- 
Mills field theory under the supposition of the spatial ho- 
mogeneity of the gauge fields [1]. For this case the La- 
grangian of pure SU (2) Yang-Mills theory reduces to 1 

L = - (A ai - ge abc A b0 A ci y - ^B 2 m {A) , (1) 

with the magnetic field B ai {A) = (l/2)ge abc e i:jk A bj A ck . 
The local SU(2) gauge invariance and the rotational in- 
variance of the original Yang-Mills action reduces to the 
symmetry under the 50(3) transformations local in time 



A^) = O(uj(t)) ab A b0 (t) - j-e abc (o(w(t))d(w(f)) 



be 



AW) = 0(u{t)) ab A bl {t) , 



(2) 



1 Everywhere in the paper we put the spatial volume V = 1. 
As result the coupling constant g becomes dimensionful with g 2 / 3 
having the dimension of energy. The volume dependence can be 
restored in the final results by replacing g 2 with g 2 /V. 



and the global spatial rotations A^ = A a jR(x)ji- 

In the constrained Hamiltonian formulation (see e.g. [2]) 
the time dependence of the gauge transformations (2) is 
exploited to put the Weyl gauge A a0 = , a=l,2,3, and 
the physical states \& have to satisfy both the Schrodinger 
equation and the three Gauss law constraints 



( 9 



G a * = -ie abc A, 



\dA ai 
8 



B 2 ai (A) 



* = E^> , (3) 



8A ci 



* = 0, a = 1,2, 3 



(4) 



The G a are the generators of the residual time indepen- 
dent gauge transformations, satisfying [G a ,H] = and 
[G a ,Gb] = ie ab cG c - Furthermore H commutes with the 
angular momentum operators Ji — —ieijkA a jd/dA a k ,i = 
1,2,3. The matrix element of an operator O is given in 
the Cartesian form 



w\om 



dA iS>'*(A)0$(A) 



(5) 



For carrying out quantum mechanical calculations it is de- 
sirable to have a corresponding unconstrained Schrodinger 
equation and to find its eigenstates in an effective way 
with high accuracy at least for the lowest states. The ba- 
sic ideas and first results of such an gauge fixed approach 
to Yang-Mills quantum mechanics will be presented in the 
following. 

2 Unconstrained Hamiltonian 

2.1 The symmetric gauge 

The local symmetry transformation (2) of the gauge po- 
tentials A a i promts us with the set of coordinates in terms 
of which the separation of the gauge degrees of freedom oc- 
curs. As in [3] I use the polar decomposition for arbitrary 
3x3 quadratic matrices 



A a i (<?, S) = Oak (q) Ski 



(G) 



with the orthogonal matrix 0(q), parametrized by the 
three angles qi, and the positive definite, symmetric 3x3 



1 



matrix S. The decomposition (6) is unique and corre- 
sponds to the symmetric gauge Xi(A) = CijkAjk=0- The 
Jacobian is \d(A ai )/d(q, S)\ oc det fl(q) Y[i<j (fa + fa), 
where fa, fa, fa are the eigenvalues of S and Qj m (q) = 
(l/2)e mH [0 T (q)dO(q)/dq j ] kl . The variables S and 
d/dS make no contribution to the Gauss law operators 
G a = —iO as (<q)Q.~ s j(q)d/dqj. Hence, assuming the invert- 
ibility of the matrix f2, the non-Abelian Gauss laws (4) 
are equivalent to the set of Abelian constraints / dqi = 
,i — 1,2,3. and the physical Hamiltonian of SU(2) 
Yang-Mills quantum mechanics reads 



H 



- & 



( 9 



\dS n 



1 



;7, 



in^rnJn + B mn (S) 



KlW-SmMl-^S))] — 



(7) 



with -fik(S) := Sik — SiktrS and the angular momenta Jj = 
—2ieijkS a jd/dS a k ,i = 1,2,3, in terms of the physical 
variables (note the factor 2). The matrix element of a 
physical operator O is given by 



(*'|0|tf) 



oc 



dS 



[Ufa 



) *'*(5)0*(5) . (8) 



The potential term in (19), has a flat valley of degener- 
ate absolute minima at zero energy, fa = fa = 0, fa — 
arbitrary, at the edge of (9). Close to the bottom of the 
valley the potential is that of a harmonic oscillator with a 
width narrowing down for larger values of fa . 

2.3 Symmetries 

As a relic of the rotational invariance of the initial gauge 
field theory the Hamiltonian (53) possesses the symmetry 



[H, J k ] = 



(13) 



with the angular momentum operators Ji = Rij£,j satis- 
fying [Ji,Jj] — itijkJk and [Ji,S,j] = . Hence the eigen- 
states of H can be characterized by the quantum numbers 
J and M. Furthermore H is invariant under arbitrary per- 
mutations Oij of any two of the three indices 1,2,3, time 
reflections T, and parity reflections P : — > — fa, 

[H,a ll } = 0, [H,T}=0, [H,P]=0. (14) 



2.4 Boundary conditions 



2.2 Unconstrained Hamiltonian in terms of ro- 
tational and scalar degrees of freedom 

In order to achieve a more transparent form for the re- 
duced Yang-Mills system (7) I shall limit myself in this 
work to the principle orbit configurations 



< fa < fa < 03 < oo 



(9) 



for the eigenvalues fa, fa, fa > of the positive definite 
symmetric matrix S (not considering singular orbits where 
two or more eigenvalues coincide) and perform a principal- 
axes transformation 

S = R(a, 0, 7 ) diag (fa, fa, fa) R T (a, (3, 7 ) , (10) 

with the SO (3) matrix R parametrized by the three Eu- 
ler angles. The Jacobian of (10) is \dS/d(a, (3, 7, fa)\ oc 
sin/3ni<j (<f>i ~ 4>j)- I n terms of the principal-axes vari- 
ables, the physical Hamiltonian reads 



The requirement of Hcrmiticity of H in the region bounded 
by the three boundary planes fa — , fa — fa , fa — fa 
and at positive infinity, leads to the conditions at fa 







V 



!>!=() 



01=0 







(15) 



corresponding to even or odd parity states respectively. I 
shall consider in the work only the first alternative allow- 
ing only even parity states in all spin sectors. Furthermore, 
due to the vanishing of the measure (12) on the boundaries 
fa — fa and fa — fa, we have to require only that 



finite A * 



finite . 



(16) 



Finally, normalisability of the wave functions requires that 
the wave functions vanish sufficiently fast at infinity 



H 



3 



1 



2 ^ 

cyclic _ 



' dfat 
Hi 



_d_ 



■77- fairr ) 2.5 Virial theorem 



+9 2 m 



(ii) 



Writing H = | (E 2 + _B 2 ) and denoting the eigenstates of 
H by \n), with energies E n , one obtains the virial theorem 



All the rotational variables in the Hamiltonian (53), which 
agrees with that of [4], are combined into the intrinsic 
angular momenta & = RfjJj obeying [&,£,•] = -itijkik- 
The matrix elements of an operator O are given as 

cyclic 



(n\E 2 \n) = 2(n\B 2 \n) , 



(17) 



() *)oc j<ln -in i,l i.h j ] [ ,lo, 



O<01<02<^3 



(12) 



It proofs to be a very useful tool to judge the quality of 
the approximate eigenstates obtained using the variational 
approach. From (17) also follows (G 2 := 2 (B 2 - E 2 )) 

(n\B 2 \n) = ^E n , (n\E 2 \n) = ^E n , (n\G 2 \n) = -^E n . 

(18) 
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3 Low energy spin-0 states 

For J 2 = £ 2 = 0, the Hamiltonian (11) reduces to 
d 2 2 / , 8 , d 



1.00- 



cyclic _ 



d 

' d<j)i 



(19) 

Similar as in the Calogero model [5] one can prove that 
the eigenstates of (19) are completely symmetric in the 
arguments (0i,0 2 ,0 3 ). We therefore can have only the 
two possible forms, the parity even 

* ( + ) (01,02,03)=* (+) (01,02,0l) , (20) 
and the parity odd 

* ( - ) (01,02,0 3 )=010203* ( - ) (0?,0|,^) , (21) 



where the functions $( ± ) are completely symmetric in the 
arguments (0 2 , 2 , </>§.)• satisfies the first and 

the second of the boundary conditions (15). I shall con- 
sider here only the even case (20), lower in energy. The 
odd case (21) can be treated completely analogously. 

3.1 Elementary symmetric polynomials 

Here it is useful to change to the new coordinates 
(ei, e2, e 3 ) defined as the elementary symmetric combina- 
tions 



ex = 0J+02 + 03 , e 2 = 0102+0203+0301 ' e 3 

The corresponding Jacobian 



9(01,02,03 



<9(ei,e 2 ,e 3 ) 



1 



(22) 



(23) 



with the square root of the discriminant 
A =]J(<P- - tf) 2 = -27el + 18 ei e 2 e 3 ~4ele 3 -4e 3 2 + e 2 2 el 

(24) 

cancels the original measure Hi<j(0i ~ 0j) m (12). Fur- 
thermore let us consider the scaling transformation 



si = ei 



S2 = 3e 2 /e 2 



S3 



27e 3 /e? , (25) 



with the Jacobian |<9(ei, e 2 , e 3 )/<9(si, s 2 , s 3 )| oc sf. Then 
the Schrodingcr equation of (19) becomes 



1 2 2 
g 5 Sl«2 



with 
where 



2si ^_- 9 A + J_/V)_!^ 
2si a s 2 V + 2 Sl ^ 4 J 

I>(°) := + D<® + D {0 1 , 



(26) 
(27) 

, (28) 
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Figure 1: The (dashed) area of integration in (s2,sg) space 
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The dotted line is the limit S3 < st/ 2 from algebraic inequali- 



ties. 



£)(0) ._ 



-12 s 2 -^-+4s 3 — + 2 
ds 2 as 3 



d 

ds 2 



(29) 
(30) 



Note that the singularities in H have disappeared via the 
transformation (22) to the elementar symmetric variables. 
The matrix elements become 



OC 

'ioi*> «y 



3/4 4 up) ( s 2) 



j 7/2 



ds 2 



ds 3 



■ 

4 up) (*2) 

ds2 J ds3 

3/4 4'™>( S2 ) 



• (31) 



Here the region of integration is given by the (s 2 , s 3 ) sat- 
isfying the inequality, obtained from 24, 

A = (-s§ + 6s 2 s 3 - 4s 3 - 4sl + 3sf) > . (32) 

Denoting the two roots of the equation A(s 3 ,S2) = 0, 
quadratic in s 3 , by 

4 up) ( S2 ) := (1 - VT^) 2 (1 + 2V1~^) , (33) 
well defined and positive for all < s 2 < 1/3, and 

4 low) ( S2 ) := (1 + VT^) 2 (1 - 2VT^T2) , (34) 

well defined for all < s 2 < 1, but positive only for 3/4 < 
s 2 < 1. The region of positive A is the shaded region in 
Fig. 1. 

3.2 Exact solution of the corresponding har- 
monic oscillator problem 

Note that for the case when the Yang-Mills potential in 
(19) is replaced by the harmonic oscillator potential 

g 2 (m + m + m) — - 2 (0? + 01. + 01) , (35) 



3 



that is g 2 S2S 2 /6 replaced by u>si/2, the above Schrodingcr 
eigenvalue problem (26) separates into a density and a de- 
formation problem 2 (k separation const., \i multiplicity), 

<S>nkn(s) = Rnk{si)Pkn(s2,S 3 ) . (36) 

The solutions of the density equations are given by 



Rnk(si) — ' 



n) 



U! 4 (U!Si) 2 



( fe -|) p -Mi/2 rfe 



_ T(n + fc+l) 
satisfying the orthonormality relations 

/>oo 

/ dsi s* /2 i2 n fe(si)i? n /fc(si) = <W , (38) 
Jo 



with the energy eigenvalues 

£„ fe = (2n + 1 + fc) w 



(39) 



The values of fc are determined by the corresponding de- 
formation problem 



(40) 



In the space of monomials s 3 , p, <? € iVo ordered 
by increasing 2p + 3g, the operator £)(°) has tridiagonal 
form, since the part is diagonal, and the parts 
and D^l act as lowering operators, as can be easily seen 
from (28)-(30). Hence can easily be diagonalised with 
eigenvectors P( p , g ) of the form 



P {p,q) = N (p,q) 



and eigenvalues 



S 2 S 3+ X] a (p^;p\q') si s\ , (41) 

2p'+3<j'<2p+3q 



fc = 2(2p + 3<z) + - . (42) 



The first 3 eigenstates P( p , q ) , ordered by increasing eigen- 
value k, are then 



k 



7/2 : P, 



(0,0) 



(43) 



fe = 15/2 : P, li0) = -J- (S2-- 



11 



221 / 2 09 / 



17 



A'2 



36 
221 



and so on. Orthonormality of the states <& with respect 
to the scalar product (31), leads to the correponding or- 
thonormality relations for the Pfc M , 

° <UP) '^ rf „„ 

(44) 



r w f 8 d *» p p 

/ ds 2 / -= Pfe' M ' Pfc M 

Jo Jo y/sE 

+ [ ds 2 [ 



i Up) ( S2 ) 



ds 3 p p 6\/3 



2 Notc that the variables si, S2 and S3 arc related to the density 
and deformation variables and7 used in [6] via si = p 2 , S2 = 
1 - /3 2 and s 3 = 1 - 3/3 2 + 2/3 3 cos 7. 



We have choosen the overall normalisation constant such 
that the constant solution, P(o,o)i corresponding to the 
lowest eigenvalue k = 7/2, is equal to one. Note that 
there is no eigenstate with k = 11/2. To the eigenvalue 
k = 31/2 correspond two states P( 3l o) an d P(o,2)- One 
can easily check that states with different values of k are 
orthonormal to each other as it should be. To orthonor- 
malise the states of equal k with respect to (45), one can 
use Gram-Schmidt orthonormalisation, P( Ptq ) — ► P^ p ^ ■ In 



(37) summary we obtain the eigenstates 

^npq(s) = Rnk(si)P[ Piq) (S2,S 3 ) 

with the energy eigenvalues 



E n ^ q = ( 2 (n + 2p + 3q) + - ) lo , 



(45) 



(46) 



which are equidistant and depend only on the total num- 
ber 2(n + 2p + 3q) of nodes. The degeneracy is therefore 
rapidly increasing with energy. 

3.3 Low energy spin-0 spectrum from varia- 
tional calculation 

To obtain the low energy spectrum of the Hamiltonian I 
shall use the Raylcigh-Ritz variational method. Lead by 
the parabolic form of the Yang-Mills potential close to the 
bottom of the classical zero-energy valley, I shall use the 
orthonormal basis <& npq in (45) as trial functions. In order 
to have rapid convergence the frequency is fixed using the 
lowest eigenstate 

$000 = \/r(9/2) w 9 / 4 exp[- WS i/2] . (47) 

The stationarity conditions for the energy functional of 
this state, 



9 9 q 2 
E[%oo] = ($oool#ol$ooo)=^ w + - 



leads to the optimal choice 

uj = \/2 g 2/3 ~ 1.259921 5 2/3 . 



(48) 



As a first upper bound for the groundstate energy of the 
Hamiltonian one therefore finds 



E < £[$000] ^ 4.25223 g 2/z . 



(49) 



In order to improve this upper bound, one truncates the 
Fock space at a certain number of nodes of the wave func- 
tions and diagonalizcs the corresponding truncated Hamil- 
tonian -fftrunk to find its eigenvalues and eigenstates. Ex- 
tending to higher and higher numbers of nodes I obtain 
the low energy spectrum in the spin-0 sector to high nu- 
merical accuracy. Including all 174(1041) trial states up 
to 30(60) nodes, I obtain for the lowest state So, 



E[S ] = 4.116719740(35) g 2/3 , 



(50) 
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(i.e. the last two digits have to be replaced by the ones in 
brackets when going from 30 to 60 nodes). The state So, 
given explicitly as (up to coefficients < 0.01) 

So = 0.98 $ooo - 0.07 $ 100 - 0.04 $ 200 

-0.16 $oio + 0.02 $ooi + 0.03 $ 02 o , (51) 

nearly coincides with the state ^ooo, The energies of the 
groundstate and the first excitations together with the cor- 
responding deviation from the virial theorem (17) 



A£ virial := 2(B 2 ) - (E 2 ) 



(52) 



are shown in Tab. 1. 

The energy spectrum of the first spin-0 eigenstates (with 
trial states up to 30 nodes) are shown in Fig. 2. The 
errorbars are inside the lines. The energies obtained here 
are in good agreement with those obtained by Luscher and 
Minister [7] and by Koller and van Baal [8], using gauge 
invariant wave functions as trial states. 

4 Higher spin states 

For non- vanishing spin we can write the Hamiltonian (11) 



i 3 



»=i 



W 1J ' k CydlC 



with the spin-0 Hamiltonian (19). The eigenstates can be 
classified according to their J 2 and J 3 quantum numbers 
and are superpositions of simultaneous eigenstates of J 2 , J3 
and £3 

\JM)= J2 *$ M ,(<t>i,<h,<h)\JMM') , (53) 

M'=0,1±,..,J± 

with the combinations (for M' > 0) 



\JMM'±) := (\JMM') ± |2M -M')) , 
v2 



(54) 



of 



\JMM') := i^^L^D%,(a,0,j) . 



Table 1: The energies and the deviation from the virial 
theorem for the lowest spin-0 eigenstates (using trial states 
up to 30 (60) nodes). 



0+ 



E[g 



273T 



4.116719740(35) 
6.386361(58) 
7.97348(34) 
9.2039(23) 
10.092(86) 
10.966(37) 
12.17(05) 



Atrial [g 2/3 ] 



W- 12 ) 

10- 10 ) 
10- 8 ) 

10- 7 ) 

10- 6 ) 

1.8- 10~ 2 (1.8 • 10- 5 ) 
-0.22 (1.5 • 10~ 4 ) 



7.6 
2.4 
4.8 
1.8 
7.5 



10~ 8 (1.2 • 
10- 7 (1.3- 

• 10- 4 (1.3 

• 10~ 3 (3.1 

• 10~ 4 (1.9 



4.1 Spin 1 

The spin-1 Schrodinger equation decays into three 
equations, one for each member of the cyclic triplet 

1*1 ' 2 ! ^3 ) — l*Ml-i *Ml+i W M0/' 



H -E+-(V 2 +V 3 ) 



= , and cycl. perm. 



One can easily show that no solutions exist which satisfy 
the boundary conditions (16). 

4.2 Spin 2 

The spin-2 Schrodinger equation decays into and into 
three equations, one for each member of the cyclic triplet 

1*1 7*2 i w 3 I — l*Ml+> W M1-' W M2-J> 



H -E+-(V 2 + V 3 ) + 2V 1 



$ 



, and cycl. perm. 



for which no solutions exist which satisfy the boundary 
conditions (16), and the coupled system 



ff -£+-(V 1 +V 2 ) 



H a -E+2V 3 + -(V 1 +V 2 ) 



(2) (2) 

for ^ M0 and ^ M2+ - The solution of this spin-2 singlet 
system (56) can be written in the form 

|2M)=tf 82, s 3 )s^Y M {(t>l, <t>\, <Pl; a, 13, 7) 

+^ 2 {s 1 ,s 2 ,s 3 )s^ 2 Y M {(t> 2 1 ,(t> 2 2 ,(t>l\a,f3, 1 ) , (56) 

with the cyclic symmetric functions ^i, 2 (s\, s 2 , S3) and 
elementary spin-2 fields Ym(</>i, 4>i, <f>3', 01, (3, 7) 



■ <p 2 )j\'2M0)+^j-(<p 1 - 2 )|2M2+) 

(57) 

, . The matrix element of 



(55) and its dual Y M := Yi 



an operator O(s) can be written as 
(2M'|0|2M)=[($' 1 |(1- S2 )0|$ 1 ) + ^(%\(s 3 - s 2 )0|$ 2 ) 

+ I(y 2 \(s3 - a 2 )0|*i> + l^' 2 \(s 2 2 - s 3 )0\* 2 ) . 

in terms of the corresponding spin-0 matrix elements (31). 
Specifying again to parity even states, the vector 
VP = ($1,^2) satisfies the Schrodinger equation of the 
same form as (26) only with the scalar replaced with 
the new matrix operators with 



D 



(2). 



/ (2(2 S2 ^+3,s 3 ^ + l)+|) 2 



\ 



(2(2^ + 3^ + 2) + !) J 
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- V ^-2-24^ J 

Again the corresponding harmonic oscillator Schrodingcr 
Equation separates 

*n/c M = i?nfc(si)P fcA1 (s2,S 3 ) . (58) 

The density equations for R n k are the same as in the spin- 
problem The values of k are determined by the corre- 
sponding deformation problem 

£ (2) P fcM (s 2 , s 3 ) = fc 2 P fcM (s 2 , 3 3 ) , (59) 

using the basis (sfsJ^O) and (0, s^a!) The first 3 eigen- 
states P(p, q ), ordered by increasing eigenvalue k, are then 

k - 11/2 : Pm-JaQ. 

= *«-i/I(T)- 

and so on. 

As for the spin-0 case one can use this orthonormal basis 
of spin-2 states of the harmonic oscillator problem as trial 
states for a variational calculation of the corresponding 
eigenstates of the spin-2 Yang-Mills Hamiltonian. The en- 
ergies of the groundstate and the first excitations together 
with the corresponding deviation from the virial theorem 
(17) are shown in Tab. 2. 

The energy spectrum of the first spin-2 eigenstates (in- 
cluding the 270 trial states with up to 30 nodes) is shown 
in Fig. 2. The errorbars are inside the lines. Also for 
spin-2 the energies obtained here are in good agreement 
with those obtained by Liischcr and Minister [7] and by 
Roller and van Baal [8]. In particular I confirm that the 
first spin-2 state is lower in energy than the first spin-0 
excitation. 



Table 2: The energies and the deviation from the virial 
theorem for the lowest spin-2 eigenstates (using trial states 
up to 30 nodes). 
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Figure 2: The energy spectrum of the first spin-0 and spin-2 
eigenstates. 
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Figure 3: The expectation values {<f>i), (4>2) and (^3) of the 
scalar fields along the principal axes as function of the energy 
for the lowest spin-0 (black boxes) and spin-2 (open circles) 
states. 

5 Expectation values of permutation 
non-invariant operators 

Expanding the obtained energy eigenstates in terms of the 
original 4>i,4>2,4>3 an d analytically calculating the Gaus- 
sian integrals in the region (9), we obtain the expectation 
values (4>i) ,(4>2) and (fe) of the three fields along the prin- 
cipal axes. For the spin-0 grounstate e.g. we obtain 

(0|^|0)+ = 0.242 g- 1 ' 3 , 
(0\<h\0)+ - 0.806 g- 1 ' 3 , 
(0|<fe|0)+ = 1.699 g- 1 ' 3 . 

Fig. 3 shows the results for the lowest spin-0 (dark sym- 
bols) and spin-2 states (open white symbols) as function of 
the energy of the states. I find that (^3) raises universally 
with increasing energy, independent whether oscillator or 
spin excitation, whereas (0i) and (^2) stay practically con- 
stant. 

The corresponding expectation values of the three in- 
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Energy [g 2/3 ] 



10 



12 



Figure 4: The expectation values (Bi), (B 2 ) and (B 3 ) of the 
magnetic fields along the principal axes as function of the en- 
ergy for the lowest spin-0 (black boxes) and spin-2 (open cir- 
cles) states. 



trinsic magnetic fields B x — g<j)2(f>3 etc. 

(0|Bi|0)+ - 1.391 5 1/3 , 
(01^10)+ = 0.415 s 1 / 3 , 
(0|B 3 |0)+ = 0.211 .g 1 / 3 , 

in accordance with the opposite ordering B\ > B 2 > B 3 , 
following from (9). The results for {Bi),{B2} and (B^) 
as function of the energy for the lowest spin-0 and spin-2 
states are shown in Fig. 4. (B 3 ) stays practically constant 
with increasing energy, whereas (-Bi) and (B^) are raising. 

Up to the considered number of 30 nodes the numerical 
errors are much smaller in both cases than the systematical 
error of the obtained states (< 1%). 

6 Summary and discussion 

It has been shown that SU(2) Yang-Mills quantum me- 
chanics can be studied very effectively and with high ac- 
curacy, by choosing the symmetric gauge and perform- 
ing a variational calculation using the orthonormal basis 
of exact eigenstates the corresponding harmonic oscillator 
problem as trial states in the spin-0 and the spin-2 sector. 
The virial theorem is shown to be fulfilled very well for 
the eigenstates obtained. The energies obtained here are 
in good agreement with those obtained by Liischer and 
Minister [7] and by Roller and van Baal [8], using gauge 
invariant wave functions as trial states. In particular I 
confirm that the first spin-2 state is lower in energy than 
the first spin-0 excitation. Furthermore I have found that 
practically all excitation energy, independently whether it 
is due to a vibrational or a rotational excitation, goes into 
the increase of the field strength along one of the three 
principal axis fields, (<ps), whereas (<f>\) and {4>2), and also 
the component (-B3) of the magnetic field along the intrin- 
sic 3-direction, remain unchanged at their vacuum values. 



Several techniques in my treatment, such as the sepa- 
ration of the harmonic oscillator problem into a density 
and a deformation problem and the diagonalization of the 
deformation operator in terms of symmetric polynomials, 
have for the case of spin-0 to some extent already been de- 
veloped in the context of the Calogero model [5] and then 
later also in the context of Yang-Mills quantum mechan- 
ics in [8]. In difference to [8], the range of integration (9) 
here is positive definite, and the corresponding range and 
measure (31) in terms of the elementary symmetric vari- 
ables is derived explicitly from the original measure (12). 
In particular the relative normalization of the deformation 
states (43) is different to that given in [8]. Furthermore 
the generalization of the use of the elementary symmetric 
variables to spin-2 and the consideration of permutation 
non-invariant expectation values here is new to the best 
of my knowledge. 
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